Research Analysis Notebook

Authors

Affiliations

Aaron Gullickson

University of Oregon

Jenifer Bratter

Rice University

This research notebook analyzes the data and models estimated for multiracial intermarriage.

Code

get_coefs <- function(x) {
  temp <- coef(x)[,c(1,2)]
  return(tibble(variable=rownames(temp), coef=temp[,1], se=temp[,2]))
}

coef_short_sym <- get_coefs(model_short_sym) |>
  filter(str_detect(variable, "race_exog")) |>
  mutate(variable=str_replace(variable, "race_exog_short_sym",""))

coef_short_asym <- get_coefs(model_short_asym) |>
  filter(str_detect(variable, "race_exog_short_asym")) |>
  mutate(variable=str_replace(variable, "race_exog_short_asym",""))

coef_full <- get_coefs(model_full_sym_multi) |>
  filter(str_detect(variable, "race_exog_full")) |>
  mutate(variable=str_replace(variable, "race_exog_full_sym",""))

coef_multi <- coef_short_sym |>
  filter(str_detect(variable, "Multi/Multi"))

Raw differences in partner selection

We begin by just looking at the raw differences in partner selection among multiracial respondents.

Code

couples <- mar_markets[[1]] |> filter(choice)

husbands <- couples |>
  mutate(race_ego=raceh, 
         race_spouse=racew,
         gender="husband") |>
  select(race_ego, race_spouse, gender)
wives <- couples |>
  mutate(race_ego=racew, 
         race_spouse=raceh,
         gender="wife") |>
  select(race_ego, race_spouse, gender)
partner_choice <- bind_rows(husbands, wives) |>
  mutate(race_spouse=str_replace(race_spouse, 
                                 "Latino, Single Race|Latino, Multiple Race",
                                 "Latino"),
         race_spouse=factor(race_spouse, 
                            levels=c("White","Black","Indigenous","Asian",
                                     "Latino","White/Black","White/Indigenous",
                                     "White/Asian","Black/Indigenous",
                                     "Black/Asian","Indigenous/Asian")))

Code

partner_choice |>
  filter(str_detect(race_ego, "/")) |>
  mutate(race_ego=factor(race_ego, 
                         levels=c("White/Asian","White/Indigenous","White/Black",
                                  "Indigenous/Asian","Black/Asian","Black/Indigenous"))) |>
  ggplot(aes(x=fct_rev(race_spouse), y=after_stat(prop), group=race_ego, 
             fill=race_ego))+
  geom_bar(position="dodge")+
  scale_y_continuous(labels=scales::percent)+
  scale_fill_manual(values = pnw_palette("Bay", 6))+
  labs(x="Race of other partner", y=NULL, fill="Race of ego partner")+
  coord_flip()+
  theme_bw()

Figure 1: Distribution of partner’s race for multiracial respondents

Unsurprisingly, for all but one multiracial group (Black/Indigenous) the most common partner race is White, which is primarily a function of group size. However, it is clear that the part-White multiracial groups have a substantially higher likelihood of being with a White partner than the non part-White multiracial groups.

By the same token, the probability of marrying another multiracial partner is low because these groups are small in the overall population. Nonetheless, we do see spikes in each of these cases for intermarriage with the exact same multiracial group, suggesting substantial endogamy among these multiracial groups.

Lets look at these same numbers in a different way by panelizing the results for each multiracial group.

Code

partner_choice |>
  filter(str_detect(race_ego, "/")) |>
  ggplot(aes(x=fct_rev(race_spouse), y=after_stat(prop), group=1))+
  geom_bar(position="dodge")+
  facet_wrap(~race_ego)+
  scale_y_continuous(labels=scales::percent)+
  scale_color_viridis_d(direction=-1, end=0.75)+
  labs(x=NULL, y=NULL)+
  coord_flip()+
  theme_bw()

Figure 2: Distribution of partner’s race for multiracial respondents, separately by panel

Two things stand out a little more clearly in Figure 2 here to me:

1. There is a tendency for each constituent group to be higher, although this can be sometimes tricky to see because of group size differences.
2. Latino is high for everybody.

Now lets do this same figure for all racial groups but separately by gender.

Code

ggplot(partner_choice,
       aes(x=fct_rev(race_spouse), y=after_stat(prop), group=gender, fill=gender))+
  geom_bar(position="dodge")+
  facet_wrap(~race_ego)+
  scale_y_continuous(labels=scales::percent)+
  scale_fill_viridis_d(direction=-1, end=0.75)+
  labs(x=NULL, y=NULL)+
  coord_flip()+
  theme_bw()

Figure 3: Distribution of partner’s race for all respondents, separately by gender and race of ego respondent

The data for part-black multiracial individuals are generally consistent with the gender pattern for black individuals - part-Black men are more likely to outmarry, particularly with Whites. For part-Asian individuals who are not also part-Black, the gender pattern seems consistent with that for Asian individuals - women are more likely than men to outmarry with Whites and less likely to marry endogamously. Its notable that the Black/Asian pattern tends to follow the Black pattern. However, its not surprising that Blackness trumps everything else.

Multiracial-to-multiracial intermarriage

Ok, now lets move to the actual model results. The first thing we want to explore is multiracial-to-multiracial intermarriage. To analyze this, we use a model that fits gender-symmetric dummies for every single multiracial-to-multiracial exogamous union. This model was difficult to fit on the full data because the size of these cells is so small, resulting in some convergence problems in the model. To resolve this issue, we created separate marriage market data only using multiracial individuals to estimate the full terms. If the models only included the racial exogamy terms, then this model would give use identical results to a model with the full data. However, because we will get slightly different estimates on the control variables in this model versus the full model, the results will differ slightly. However, they should be sufficiently close to allow us to determine a more reduced form multiracial-to-multiracial coding we can use in the full data models.

Figure 4 below shows the odds ratios of outmarriage to each other multiracial group for each ego multiracial group, separately by panel.

Code

# the somewhat odd regular expression \\..*/ will identify all the names with a 
# period followed by some characters followed by a forward slash. This will give
# me all cases of a single race/multiple race marriage
temp <- coef_full |>
  filter(str_detect(variable, "/.*\\.*/")) |>
  mutate(race_ego=str_split(variable, "\\.", simplify=TRUE)[,1],
         race_choice=str_split(variable, "\\.", simplify=TRUE)[,2],
         race_cons1=str_split(race_ego, "/", simplify=TRUE)[,1],
         race_cons2=str_split(race_ego, "/", simplify=TRUE)[,2],
         overlap=str_detect(race_choice, race_cons1) | 
           str_detect(race_choice, race_cons2)) |>
  select(race_ego, race_choice, overlap, coef, se)
         
# this needs to be duplicated so that each possibility is both ego and choice.
# easiest way to do this is to create a second one, flip the variables and 
# bind them together
temp2 <- temp
temp <- temp2 |>
  mutate(race_ego=temp$race_choice,
         race_choice=temp$race_ego) |>
  bind_rows(temp) |>
  mutate(race_ego=factor(race_ego, 
                         levels=c("White/Asian","Black/Asian",
                                  "White/Indigenous","Black/Indigenous",
                                  "White/Black","Indigenous/Asian")),
         race_choice=factor(race_choice, 
                         levels=c("White/Asian","Black/Asian",
                                  "White/Indigenous","Black/Indigenous",
                                  "White/Black","Indigenous/Asian")))


ggplot(temp, aes(x=race_choice, y=exp(coef), color=overlap))+
  geom_linerange(linewidth=1.25, 
                 aes(ymax=exp(coef+1.386*se), ymin=exp(coef-1.386*se)))+
  geom_linerange(linewidth=0.5, 
                 aes(ymax=exp(coef+1.96*se), ymin=exp(coef-1.96*se)))+
  geom_point(size=2)+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  coord_flip()+
  facet_wrap(~race_ego, ncol=2)+
  scale_color_viridis_d(direction=-1, end=0.75)+
  labs(x="multiple race partner", 
       y="odds of marriage with multiple race partner relative to multiracial endogamy",
       color="partial noverlap")+
  theme_bw()

Figure 4: Odds of interracial marriage with other multiracial groups for each multiracial group. The purple dots indicate cases where the ego multiracial group has overlapping partial ancestry with the other multiracial group, while the green dots indicate cases where is no overlap. 83.4% and 95% confidence intervals are shown for each estimate with thick and thin bars, respectively. The grey line indicates the case of an odds equal to one, meaning this type of intermarriage is as likely as endogamy (e.g. White/Black person marrying White/Black person).

As expected, some of the confidence intervals are pretty large because we have smaller samples of these types of marriages. I notice two things:

1. None of these odds ratios are particularly high. The point estimates never break the 0.25 level. Below, we will calculate some of these same odds for multiracial-to-monoracial and for monoracial-to-monoracial. These odds are consistent with some of the higher monoracial-to-monoracial odds ratios but not outstanding. Thus, there is no evidence of breaking through to some sort of post-racial multiracial eve utopia here. Race still matters even if in combination.
2. The results tend to show that when the two multiracial have a partial overlap in ancestry (e.g. White/Asian and White/Black share White), as indicated by the coloring, the odds tend to be higher.

To get another view of the results, lets line them all up from highest to lowest on the same plot.

Code

coef_full |>
  mutate(race_ego=str_split(variable, "\\.", simplify=TRUE)[,1],
         race_choice=str_split(variable, "\\.", simplify=TRUE)[,2],
         race_cons1=str_split(race_ego, "/", simplify=TRUE)[,1],
         race_cons2=str_split(race_ego, "/", simplify=TRUE)[,2],
         overlap=str_detect(race_choice, race_cons1) | 
           str_detect(race_choice, race_cons2),
         variable=str_replace(variable, "\\.","+")) |>
  ggplot(aes(x=reorder(variable, coef, max), y=exp(coef), color=overlap))+
  geom_linerange(linewidth=1.25, 
                 aes(ymax=exp(coef+1.386*se), ymin=exp(coef-1.386*se)))+
  geom_linerange(linewidth=0.5, 
                 aes(ymax=exp(coef+1.96*se), ymin=exp(coef-1.96*se)))+
  geom_point(size=2)+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  coord_flip()+
  scale_color_viridis_d(direction=-1, end=0.75)+
  labs(x=NULL, 
       y="odds of exogamous marriage with multiple race partner relative to multiracial endogamy",
       color="partial overlap")+
  theme_bw()

Figure 5: Odds of interracial marriage with other multiracial groups for each multiracial group, ordered from highest to lowest. The purple dots indicate cases where the ego multiracial group has overlapping partial ancestry with the other multiracial group, while the green dots indicate cases where is no overlap. 83.4% and 95% confidence intervals are shown for each estimate with thick and thin bars, respectively. The grey line indicates the case of an odds equal to one, meaning this type of intermarriage is as likely as endogamy (e.g. White/Black person marrying White/Black person). The dashed colored lines indicate the results from a model with two parameters for shared and not shared ancestry.

The three non-overlapping cases are at the bottom, although Black/Indigenous+Indigenous/Asian is also there (with huge error bars). A simple way to reduce the terms for this part of the model would just be to separate non-overlapping from overlapping cases, which would reduce us from 15 terms to 2 terms. Although there is some substantial variation in addition to this, I cannot see a simple pattern that could account for it, especially with the error bars.

Multiracial-to-monoracial intermarriage

Figure 6 shows the odds of intermarriage between each multiracial group and each single race group. For reference, I also draw lines indicating the reduced form odds of intermarriage between multiracial groups (shared vs. no shared ancestry).

Code

# the somewhat odd regular expression \\..*/ will identify all the names with a 
# period followed by some characters followed by a forward slash. This will give
# me all cases of a single race/multiple race marriage
coef_short_sym |>
  filter(str_detect(variable, "\\..*/")) |>
  mutate(single_race=str_split(variable, "\\.", simplify=TRUE)[,1],
         multiple_race=str_split(variable, "\\.", simplify=TRUE)[,2],
         constituent=str_detect(multiple_race, single_race),
         single_race=factor(single_race, 
                            levels=c("Latino, Multiple Race", 
                                     "Latino, Single Race","Asian","Indigenous",
                                     "Black","White")),
         multiple_race=factor(multiple_race,
                              levels=c("White/Asian","Black/Asian",
                                       "White/Indigenous","Black/Indigenous",
                                       "White/Black","Indigenous/Asian")),
         upper84=exp(coef+1.386*se),
         upper95=exp(coef+1.96*se),
         lower84=exp(coef-1.386*se),
         lower95=exp(coef-1.96*se)) |>
  ggplot(aes(x=single_race, y=exp(coef), color=constituent))+
  geom_hline(data=coef_multi, aes(yintercept=exp(coef), linetype=variable),
             color="red")+
  geom_linerange(linewidth=1.25, 
                  aes(ymax=upper84, ymin=lower84))+
  geom_linerange(linewidth=0.5, 
                  aes(ymax=upper95, ymin=lower95))+
  geom_point(size=2)+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  coord_flip()+
  facet_wrap(~multiple_race, ncol=2)+
  scale_color_viridis_d(direction=-1, end=0.75)+
  labs(x="single race partner", 
       y="odds of marriage with single race partner relative to multiracial endogamy")+
  theme_bw()

Figure 6: Odds of interracial marriage with single race groups for each multiracial group. The purple dots indicate cases where the single race group is constituent of the multiracial group, while the green dots indicate cases where the single race group is not constituent of the multiracial group. 83.4% and 95% confidence intervals are shown for each estimate with thick and thin bars, respectively. The orange dashed line indicates the odds of intermarriage to a multiracial individual with no shared ancestry and the blue dashed line indicates the odds of intermarriage to a multiracial individual with shared ancestry. The grey line indicates the case of an odds equal to one, meaning this type of intermarriage is as likely as endogamy.

A few things I note about Figure 6:

1. Latinos have a high rate of intermarriage with all the multiracial groups. This is not accounted for simply by multiple vs. single race Latinos.
2. Notwithstanding (1), the constituent cases tend of have higher odds. This is most obvious in the case of White/Asian and White/Black partners. The two deviations from this are:

• White/Indigenous - the odds of intermarriage with an Indigenous person are as low as the odds for Blacks and Asians.
• Indigenous/Asians - again the odds of intermarriage with an Indigenous person are relatively low and the odds of intermarriage with a White person is the highest. However, this group is also the noisiest because it is the smallest.

3. For most groups the odds of multiracial-to-multiracial intermarriage with shared ancestry are close to the odds for the highest single race group, but not substantially higher. For White/Asians, they are substantially lower.

fig-odds-ratios-latino takes a closer look at the results for Latinos.

Code

coef_short_sym |>
  filter(str_detect(variable, "Latino")) |>
  mutate(multiple=case_when(
    str_detect(variable, "Multiple Race") ~ "Multiple race",
    TRUE ~ "Single race"),
    str_detect(variable, "Multiple Race"),
    variable=str_replace(variable, "Latino,\\s.+\\sRace",""),
    variable=str_replace(variable, "\\.", ""),
    type=case_when(
      str_count(variable, "/")>0 ~ "Multiracial partner",
      TRUE ~ "Monoracial partner"
    )) |>
  filter(variable!="") |>
  ggplot(aes(x=reorder(variable, coef, max), y=exp(coef), 
             shape=multiple, color=type))+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  geom_point()+
  scale_color_viridis_d(direction=-1, end=0.75)+
  labs(x="Race of non-Latino partner",
       y="odds of intermarriage for Latino partners relative to endogamy",
    color="non-latino partner", shape="Latino partner")+
  coord_flip()+
  theme_bw()

Figure 7: Odds of interracial marriage with single and multiracial race groups for Latinos Results are shown separately for Latinos who selected a single race and those who selected multiple races. The grey line indicates the case of an odds equal to one, meaning this type of intermarriage is as likely as endogamy.

In general, multiracial Latino partner were slightly more likely to intermarry with all of the multiracial groups (except Black/Indigenous) than single race Latino partners, and vice versa for single race groups. However, these differences were not large in most cases. The exception is intermarriage with Whites which is far less common among multiple race Latinos. Aside from intermarriage with whites, Latino intermarriage is more likely to non-Latino multiracial groups than non-Latino single race groups.

To get at potentially gender differences, Figure 8 displays a similar graph but with separate results depending on the gender combination from a model that allows for gender asymmetry.

Code

# before period is husband, after period is wife
temp <- coef_short_asym |>
  filter(str_detect(variable, "/") & 
           !str_detect(variable, "Multi/Multi|Latino")) |>
  mutate(race_husband=str_split(variable, "\\.", simplify=TRUE)[,1],
         race_wife=str_split(variable, "\\.", simplify=TRUE)[,2],
         race_label=ifelse(str_detect(race_husband, "/"), 
                           paste(race_wife, " Female", 
                                 sep=""),
                           paste(race_husband, " Male", 
                                 sep="")),
         multiple_race=ifelse(str_detect(race_husband, "/"), race_husband,
                            race_wife),
         single_race=ifelse(str_detect(race_husband, "/"), race_wife,
                            race_husband),
         race_group=paste(single_race, ".", multiple_race, sep=""),
         single_race=factor(single_race, 
                            levels=c("Latino, Multiple Race", 
                                     "Latino, Single Race","Asian","Indigenous",
                                     "Black","White")),
         multiple_race=factor(multiple_race,
                              levels=c("White/Asian","Black/Asian",
                                       "White/Indigenous","Black/Indigenous",
                                       "White/Black","Indigenous/Asian")),
         upper84=exp(coef+1.386*se),
         upper95=exp(coef+1.96*se),
         lower84=exp(coef-1.386*se),
         lower95=exp(coef-1.96*se)) |>
  arrange(multiple_race, single_race, race_label)

temp |>
  #filter(race_group!="Indigenous.Black/Indigenous" &  # CIs too big
  #         race_group!="Indigenous.Black/Asian") |>
  ggplot(aes(x=fct_inorder(race_label), y=exp(coef), color=single_race))+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  geom_linerange(linewidth=1.25, 
                 aes(ymax=upper84, ymin=lower84))+
  geom_linerange(linewidth=0.5, 
                 aes(ymax=upper95, ymin=lower95))+
  geom_point(size=2)+
  coord_flip()+
  facet_wrap(~multiple_race, ncol=2)+
  scale_color_viridis_d(direction=-1, end=0.9)+
  theme_bw()+
  labs(color="Single race partner",
       x="single race partner", 
       y="odds of marriage with single race partner relative to multiracial endogamy")

Figure 8: Odds of interracial marriage with single race groups for each multiracial group, separately by the gender combination of the partners. 83.4% and 95% confidence intervals are shown for each estimate with thick and thin bars, respectively. The grey line indicates the case of an odds equal to one, meaning this type of intermarriage is as likely as endogamy. Some upper confidence intervals are truncated to maintain focus on point estimates.

What do we note from Figure 8?

1. For all part-Black individuals, the pattern is like that for single race Black individuals. part-Black male/non-Black female pairings are more likely than vice versa, and correspondingly unions between part-Black female/Black male are more common than vice versa. This is easiest to see for the White/Black case. This implies that with regard to gender dynamics, the one-drop rule operates for people of part-Black ancestry.
2. We don’t really see the same for people of part-Asian ancestry. It is marginally more likely to see White/Asian female+White male pairings than vice versa but the difference is not statistically significant. Furthermore, White/Asians have the same gendered pattern of intermarriage with Asians as other groups - women less likely to marry an Asian male and more likely to marry an Asian female. This suggests more flexibility in how partially Asian people are treated on the marriage market.

Monoracial-to-monoracial intermarriage

We don’t care as much about the monoracial-to-monoracial comparisons, but lets go ahead and plot those out since we have them.

Code

coef_short_sym |>
  filter(!str_detect(variable, "/")) |> 
  mutate(variable=str_replace(variable, "\\.","+"),
         upper84=exp(coef+1.386*se),
         upper95=exp(coef+1.96*se),
         lower84=exp(coef-1.386*se),
         lower95=exp(coef-1.96*se)) |>
  ggplot(aes(x=reorder(variable, coef, max), y=exp(coef)))+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  geom_linerange(linewidth=1.25, color="grey20",
                 aes(ymax=upper84, ymin=lower84))+
  geom_linerange(linewidth=0.5, color="grey20",
                 aes(ymax=upper95, ymin=lower95))+
  geom_point(color="grey20")+
  coord_flip()+
  theme_bw()+
  labs(x=NULL, 
       y="odds of interracial union relative to endogamy")

Figure 9: Odds of interracial marriage between single race groups. 83.4% and 95% confidence intervals are shown for each estimate with thick and thin bars, respectively. The grey line indicates the case of an odds equal to one, meaning this type of intermarriage is as likely as endogamy.

A few things worth nothing:

1. All of the least likely cases involve a Black partner, except for Black+Latino intermarriage.
2. White/Latino is far the most likely combo, but only for Latino, Single Race.
3. In general, Latino multiple race partners are less exogamous than Latino single race respondents.

For a better comparison across all cases, lets do a megaplot to get everything on the board at once:

Code

coef_short_sym |>
  mutate(variable=str_replace(variable, "\\.","+"),
         upper84=exp(coef+1.386*se),
         upper95=exp(coef+1.96*se),
         lower84=exp(coef-1.386*se),
         lower95=exp(coef-1.96*se),
         type=case_when(
           str_detect(variable, "Multi/Multi") ~ "Multi-to-Multi",
           str_count(variable, "/")>0 ~ "Multi-to-Mono",
           TRUE ~ "Mono-to-Mono"
         )) |>
  ggplot(aes(x=reorder(variable, coef, max), y=exp(coef), color=type))+
  geom_hline(yintercept=1, linewidth=0.5, color="grey30")+
  geom_linerange(linewidth=1.25,
                 aes(ymax=upper84, ymin=lower84))+
  geom_linerange(linewidth=0.5,
                 aes(ymax=upper95, ymin=lower95))+
  geom_point()+
  coord_flip()+
  scale_color_viridis_d(direction=-1, end=0.75)+
  theme_bw()+
  theme(legend.position = "bottom")+
  labs(x=NULL, 
       y="odds of interracial union relative to endogamy")

Figure 10: Later

Full Models

Code

htmlreg(lapply(list(model_full_sym_multi, model_short_sym_multi,
                      model_short_sym, model_short_asym), convertModel),
        digits=3,
        caption="Model predicting union formation",
        caption.above=TRUE)

?(caption)

Model predicting union formation

	Model 1	Model 2	Model 3	Model 4
agediff	0.063***	0.062***	0.071***	0.070***
	(0.006)	(0.006)	(0.000)	(0.000)
agediff^2	-0.011***	-0.011***	-0.012***	-0.012***
	(0.001)	(0.001)	(0.000)	(0.000)
hypergamyTRUE	-0.080	-0.079	-0.047***	-0.036***
	(0.220)	(0.209)	(0.007)	(0.007)
hypogamyTRUE	0.258	0.261	0.243***	0.237***
	(0.157)	(0.151)	(0.008)	(0.008)
edcross_hsTRUE	-0.653***	-0.653***	-0.892***	-0.897***
	(0.162)	(0.161)	(0.008)	(0.008)
edcross_scTRUE	-0.565***	-0.564***	-0.674***	-0.676***
	(0.137)	(0.136)	(0.006)	(0.006)
edcross_cTRUE	-0.765***	-0.752***	-0.821***	-0.824***
	(0.141)	(0.136)	(0.006)	(0.006)
birth_endogTRUE	0.344	0.352	0.406***	0.410***
	(0.216)	(0.227)	(0.015)	(0.016)
language_endogTRUE	1.603***	1.576***	1.414***	1.419***
	(0.241)	(0.251)	(0.008)	(0.008)
race_altWhite/Indigenous	0.234	0.293	0.958***	0.961***
	(0.174)	(0.172)	(0.052)	(0.054)
race_altWhite/Asian	0.012	0.100	0.852***	0.854***
	(0.161)	(0.149)	(0.113)	(0.116)
race_altBlack/Indigenous	1.014*	0.579*	1.467***	1.473***
	(0.408)	(0.245)	(0.268)	(0.266)
race_altBlack/Asian	0.768*	0.781*	1.746***	1.757***
	(0.344)	(0.359)	(0.259)	(0.271)
race_altIndigenous/Asian	0.548*	0.584**	1.043***	1.049***
	(0.246)	(0.220)	(0.208)	(0.208)
race_exog_full_symWhite/Black.White/Indigenous	-2.148***
	(0.174)
race_exog_full_symWhite/Black.White/Asian	-1.854***
	(0.177)
race_exog_full_symWhite/Black.Black/Indigenous	-2.747***
	(0.350)
race_exog_full_symWhite/Black.Black/Asian	-1.467**
	(0.480)
race_exog_full_symWhite/Black.Indigenous/Asian	-3.033***
	(0.768)
race_exog_full_symWhite/Indigenous.White/Asian	-1.440***
	(0.123)
race_exog_full_symWhite/Indigenous.Black/Indigenous	-2.677***
	(0.358)
race_exog_full_symWhite/Indigenous.Black/Asian	-3.094***
	(0.574)
race_exog_full_symWhite/Indigenous.Indigenous/Asian	-1.526***
	(0.286)
race_exog_full_symWhite/Asian.Black/Indigenous	-3.999***
	(0.526)
race_exog_full_symWhite/Asian.Black/Asian	-2.148***
	(0.473)
race_exog_full_symWhite/Asian.Indigenous/Asian	-2.166***
	(0.350)
race_exog_full_symBlack/Indigenous.Black/Asian	-2.353***
	(0.589)
race_exog_full_symBlack/Indigenous.Indigenous/Asian	-3.074*
	(1.384)
race_exog_full_symBlack/Asian.Indigenous/Asian	-2.394*
	(1.152)
race_exog_short_symMulti/Multi, no shared ancestry		-3.266***	-3.154***
		(0.380)	(0.370)
race_exog_short_symMulti/Multi, shared ancestry		-1.871***	-1.695***
		(0.091)	(0.142)
race_altBlack			0.249***	0.245***
			(0.011)	(0.011)
race_altIndigenous			0.941***	0.943***
			(0.032)	(0.032)
race_altAsian			0.917***	0.931***
			(0.026)	(0.026)
race_altLatino, Single Race			0.455***	0.456***
			(0.010)	(0.010)
race_altLatino, Multiple Race			1.411***	1.412***
			(0.056)	(0.056)
race_altWhite/Black			0.825***	0.825***
			(0.107)	(0.109)
race_exog_short_symAsian.Black/Asian			-2.312***
			(0.320)
race_exog_short_symAsian.Black/Indigenous			-4.064***
			(0.470)
race_exog_short_symAsian.Indigenous/Asian			-1.784***
			(0.334)
race_exog_short_symAsian.Latino, Multiple Race			-2.778***
			(0.114)
race_exog_short_symAsian.Latino, Single Race			-2.269***
			(0.042)
race_exog_short_symAsian.White/Asian			-1.389***
			(0.157)
race_exog_short_symAsian.White/Black			-2.635***
			(0.171)
race_exog_short_symAsian.White/Indigenous			-2.698***
			(0.140)
race_exog_short_symBlack.Asian			-3.253***
			(0.055)
race_exog_short_symBlack.Black/Asian			-1.884***
			(0.285)
race_exog_short_symBlack.Black/Indigenous			-1.915***
			(0.257)
race_exog_short_symBlack.Indigenous			-3.393***
			(0.078)
race_exog_short_symBlack.Indigenous/Asian			-3.075***
			(0.337)
race_exog_short_symBlack.Latino, Multiple Race			-2.807***
			(0.086)
race_exog_short_symBlack.Latino, Single Race			-2.303***
			(0.026)
race_exog_short_symBlack.White/Asian			-3.276***
			(0.157)
race_exog_short_symBlack.White/Black			-1.468***
			(0.120)
race_exog_short_symBlack.White/Indigenous			-3.789***
			(0.117)
race_exog_short_symIndigenous.Asian			-3.012***
			(0.147)
race_exog_short_symIndigenous.Black/Asian			-4.295***
			(0.749)
race_exog_short_symIndigenous.Black/Indigenous			-3.045***
			(0.413)
race_exog_short_symIndigenous.Indigenous/Asian			-2.630***
			(0.373)
race_exog_short_symIndigenous.Latino, Multiple Race			-3.289***
			(0.213)
race_exog_short_symIndigenous.Latino, Single Race			-2.210***
			(0.066)
race_exog_short_symIndigenous.White/Asian			-3.269***
			(0.212)
race_exog_short_symIndigenous.White/Black			-3.322***
			(0.255)
race_exog_short_symIndigenous.White/Indigenous			-2.744***
			(0.127)
race_exog_short_symLatino, Multiple Race.Black/Asian			-2.402***
			(0.546)
race_exog_short_symLatino, Multiple Race.Black/Indigenous			-2.243***
			(0.554)
race_exog_short_symLatino, Multiple Race.Indigenous/Asian			-2.304***
			(0.355)
race_exog_short_symLatino, Multiple Race.White/Asian			-1.698***
			(0.216)
race_exog_short_symLatino, Multiple Race.White/Black			-1.449***
			(0.243)
race_exog_short_symLatino, Multiple Race.White/Indigenous			-1.631***
			(0.177)
race_exog_short_symLatino, Single Race.Black/Asian			-2.598***
			(0.303)
race_exog_short_symLatino, Single Race.Black/Indigenous			-2.932***
			(0.302)
race_exog_short_symLatino, Single Race.Indigenous/Asian			-2.040***
			(0.283)
race_exog_short_symLatino, Single Race.Latino, Multiple Race			-2.082***
			(0.073)
race_exog_short_symLatino, Single Race.White/Asian			-1.653***
			(0.088)
race_exog_short_symLatino, Single Race.White/Black			-1.704***
			(0.128)
race_exog_short_symLatino, Single Race.White/Indigenous			-1.876***
			(0.101)
race_exog_short_symWhite.Asian			-1.627***
			(0.021)
race_exog_short_symWhite.Black			-3.102***
			(0.014)
race_exog_short_symWhite.Black/Asian			-3.034***
			(0.278)
race_exog_short_symWhite.Black/Indigenous			-3.534***
			(0.307)
race_exog_short_symWhite.Indigenous			-1.993***
			(0.036)
race_exog_short_symWhite.Indigenous/Asian			-1.730***
			(0.239)
race_exog_short_symWhite.Latino, Multiple Race			-2.179***
			(0.065)
race_exog_short_symWhite.Latino, Single Race			-1.126***
			(0.010)
race_exog_short_symWhite.White/Asian			-1.163***
			(0.102)
race_exog_short_symWhite.White/Black			-1.778***
			(0.107)
race_exog_short_symWhite.White/Indigenous			-1.501***
			(0.057)
race_exog_short_asymAsian.Black				-4.216***
				(0.132)
race_exog_short_asymAsian.Black/Asian				-2.947***
				(0.427)
race_exog_short_asymAsian.Black/Indigenous				-5.547***
				(1.093)
race_exog_short_asymAsian.Indigenous				-3.208***
				(0.162)
race_exog_short_asymAsian.Indigenous/Asian				-1.969***
				(0.382)
race_exog_short_asymAsian.Latino, Multiple Race				-3.244***
				(0.183)
race_exog_short_asymAsian.Latino, Single Race				-2.497***
				(0.050)
race_exog_short_asymAsian.White				-2.168***
				(0.025)
race_exog_short_asymAsian.White/Asian				-1.783***
				(0.176)
race_exog_short_asymAsian.White/Black				-3.627***
				(0.265)
race_exog_short_asymAsian.White/Indigenous				-3.298***
				(0.211)
race_exog_short_asymBlack.Asian				-2.746***
				(0.079)
race_exog_short_asymBlack.Black/Asian				-1.583***
				(0.334)
race_exog_short_asymBlack.Black/Indigenous				-1.633***
				(0.274)
race_exog_short_asymBlack.Indigenous				-2.939***
				(0.088)
race_exog_short_asymBlack.Indigenous/Asian				-2.510***
				(0.353)
race_exog_short_asymBlack.Latino, Multiple Race				-2.396***
				(0.097)
race_exog_short_asymBlack.Latino, Single Race				-1.911***
				(0.040)
race_exog_short_asymBlack.White				-2.669***
				(0.015)
race_exog_short_asymBlack.White/Asian				-2.909***
				(0.164)
race_exog_short_asymBlack.White/Black				-1.228***
				(0.134)
race_exog_short_asymBlack.White/Indigenous				-3.334***
				(0.140)
race_exog_short_asymBlack/Asian.Asian				-1.914***
				(0.383)
race_exog_short_asymBlack/Asian.Black				-2.342***
				(0.385)
race_exog_short_asymBlack/Asian.Indigenous				-3.788***
				(0.957)
race_exog_short_asymBlack/Asian.Latino, Multiple Race				-2.256***
				(0.616)
race_exog_short_asymBlack/Asian.Latino, Single Race				-2.468***
				(0.317)
race_exog_short_asymBlack/Asian.White				-2.931***
				(0.286)
race_exog_short_asymBlack/Indigenous.Asian				-3.491***
				(0.458)
race_exog_short_asymBlack/Indigenous.Black				-2.284***
				(0.296)
race_exog_short_asymBlack/Indigenous.Indigenous				-2.631***
				(0.431)
race_exog_short_asymBlack/Indigenous.Latino, Multiple Race				-2.219***
				(0.604)
race_exog_short_asymBlack/Indigenous.Latino, Single Race				-2.687***
				(0.305)
race_exog_short_asymBlack/Indigenous.White				-3.407***
				(0.286)
race_exog_short_asymIndigenous.Asian				-2.859***
				(0.198)
race_exog_short_asymIndigenous.Black				-4.122***
				(0.128)
race_exog_short_asymIndigenous.Black/Asian				-5.056**
				(1.574)
race_exog_short_asymIndigenous.Black/Indigenous				-3.652***
				(0.626)
race_exog_short_asymIndigenous.Indigenous/Asian				-2.446***
				(0.464)
race_exog_short_asymIndigenous.Latino, Multiple Race				-3.324***
				(0.292)
race_exog_short_asymIndigenous.Latino, Single Race				-2.312***
				(0.102)
race_exog_short_asymIndigenous.White				-1.967***
				(0.049)
race_exog_short_asymIndigenous.White/Asian				-3.310***
				(0.311)
race_exog_short_asymIndigenous.White/Black				-3.639***
				(0.361)
race_exog_short_asymIndigenous.White/Indigenous				-2.552***
				(0.151)
race_exog_short_asymIndigenous/Asian.Asian				-1.603***
				(0.374)
race_exog_short_asymIndigenous/Asian.Black				-4.302***
				(0.749)
race_exog_short_asymIndigenous/Asian.Indigenous				-2.733***
				(0.741)
race_exog_short_asymIndigenous/Asian.Latino, Multiple Race				-2.189***
				(0.560)
race_exog_short_asymIndigenous/Asian.Latino, Single Race				-2.145***
				(0.285)
race_exog_short_asymIndigenous/Asian.White				-1.920***
				(0.232)
race_exog_short_asymLatino, Multiple Race.Asian				-2.417***
				(0.113)
race_exog_short_asymLatino, Multiple Race.Black				-3.469***
				(0.130)
race_exog_short_asymLatino, Multiple Race.Black/Asian				-2.711**
				(0.830)
race_exog_short_asymLatino, Multiple Race.Black/Indigenous				-2.184*
				(1.003)
race_exog_short_asymLatino, Multiple Race.Indigenous				-3.236***
				(0.327)
race_exog_short_asymLatino, Multiple Race.Indigenous/Asian				-2.388***
				(0.485)
race_exog_short_asymLatino, Multiple Race.Latino, Single Race				-2.063***
				(0.080)
race_exog_short_asymLatino, Multiple Race.White				-2.152***
				(0.085)
race_exog_short_asymLatino, Multiple Race.White/Asian				-1.602***
				(0.271)
race_exog_short_asymLatino, Multiple Race.White/Black				-1.556***
				(0.294)
race_exog_short_asymLatino, Multiple Race.White/Indigenous				-1.292***
				(0.227)
race_exog_short_asymLatino, Single Race.Asian				-2.086***
				(0.052)
race_exog_short_asymLatino, Single Race.Black				-2.901***
				(0.038)
race_exog_short_asymLatino, Single Race.Black/Asian				-2.759***
				(0.393)
race_exog_short_asymLatino, Single Race.Black/Indigenous				-3.309***
				(0.386)
race_exog_short_asymLatino, Single Race.Indigenous				-2.120***
				(0.065)
race_exog_short_asymLatino, Single Race.Indigenous/Asian				-1.952***
				(0.378)
race_exog_short_asymLatino, Single Race.Latino, Multiple Race				-2.105***
				(0.083)
race_exog_short_asymLatino, Single Race.White				-1.121***
				(0.012)
race_exog_short_asymLatino, Single Race.White/Asian				-1.754***
				(0.104)
race_exog_short_asymLatino, Single Race.White/Black				-1.970***
				(0.137)
race_exog_short_asymLatino, Single Race.White/Indigenous				-1.853***
				(0.114)
race_exog_short_asymMulti/Multi, no shared ancestry				-3.155***
				(0.371)
race_exog_short_asymMulti/Multi, shared ancestry				-1.698***
				(0.144)
race_exog_short_asymWhite.Asian				-1.293***
				(0.027)
race_exog_short_asymWhite.Black				-3.739***
				(0.023)
race_exog_short_asymWhite.Black/Asian				-3.166***
				(0.345)
race_exog_short_asymWhite.Black/Indigenous				-3.727***
				(0.379)
race_exog_short_asymWhite.Indigenous				-2.020***
				(0.037)
race_exog_short_asymWhite.Indigenous/Asian				-1.583***
				(0.279)
race_exog_short_asymWhite.Latino, Multiple Race				-2.207***
				(0.064)
race_exog_short_asymWhite.Latino, Single Race				-1.129***
				(0.012)
race_exog_short_asymWhite.White/Asian				-1.112***
				(0.090)
race_exog_short_asymWhite.White/Black				-1.990***
				(0.122)
race_exog_short_asymWhite.White/Indigenous				-1.459***
				(0.067)
race_exog_short_asymWhite/Asian.Asian				-1.088***
				(0.163)
race_exog_short_asymWhite/Asian.Black				-3.836***
				(0.231)
race_exog_short_asymWhite/Asian.Indigenous				-3.206***
				(0.345)
race_exog_short_asymWhite/Asian.Latino, Multiple Race				-1.807***
				(0.237)
race_exog_short_asymWhite/Asian.Latino, Single Race				-1.560***
				(0.104)
race_exog_short_asymWhite/Asian.White				-1.225***
				(0.126)
race_exog_short_asymWhite/Black.Asian				-2.019***
				(0.220)
race_exog_short_asymWhite/Black.Black				-1.875***
				(0.124)
race_exog_short_asymWhite/Black.Indigenous				-2.983***
				(0.344)
race_exog_short_asymWhite/Black.Latino, Multiple Race				-1.353***
				(0.283)
race_exog_short_asymWhite/Black.Latino, Single Race				-1.443***
				(0.146)
race_exog_short_asymWhite/Black.White				-1.561***
				(0.105)
race_exog_short_asymWhite/Indigenous.Asian				-2.335***
				(0.156)
race_exog_short_asymWhite/Indigenous.Black				-4.481***
				(0.194)
race_exog_short_asymWhite/Indigenous.Indigenous				-2.952***
				(0.209)
race_exog_short_asymWhite/Indigenous.Latino, Multiple Race				-1.987***
				(0.210)
race_exog_short_asymWhite/Indigenous.Latino, Single Race				-1.896***
				(0.113)
race_exog_short_asymWhite/Indigenous.White				-1.546***
				(0.060)
Deviance	4019	4060	1579219	1575065
BIC (relative to null)	-3909	-3966	-1869518	-1872976
***p < 0.001; **p < 0.01; *p < 0.05